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Self-consistency of interaction models QGSJET 01, SIBYLL 2.1, NEXUS 3.97 and QGSJET II is checked in 
terms of their ability to reproduce simultaneously experimental data on fluxes of muons and hadrons. From this 
point of view SIBYLL 2.1 gives the most acceptable, though not quite satisfactory, results. Analysis of the situation 
for muons supports our previous conclusions, that high-energy muon deficit is due both to underestimation of 
primary light nuclei fluxes in direct emulsion chamber experiments and to softness of p + A — > , K ± + X 
inclusive spectra in fragmentation region, especially prominent in case of QGSJET 01 model. 



1. Introduction 

At present information on the characteristics 
of hadronic interactions in fragmentation region 
is still scarce or missing and experiments with 
'roman pots' are anticipated to improve the sit- 
uation. Some of this information, in principle, 
could be obtained with the use of the data on CR 
muon and hadron spectra, provided PCR spectra 
are known with high precision, but that is not the 
case. The obvious obstacle here is that at high 
energies PCR fluxes themselves are functionals 
of various interaction parameters plus their accu- 
racy is appreciably affected by additional system- 
atic effects. However, comparison of the hadron 
and muon fluxes, predicted by different interac- 
tion packages with the experimental data still al- 
lows to get information on the fragmentation par- 
ticle spectra in the quasi-independent on the PCR 
fluxes data way. Besides, as shown below in this 
paper, our notions on behavior of the PCR light 
nuclei spectra can also benefit from such analysis. 
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Figure 1. Muon spectra at sea level for PCR fits 
from [1]. Experimental data: L3+C [2], LVD [3], 
underground — the lower bound of underground 
measurements [4]. 



2. Muons 

Average numbers of hadrons and muons in EAS 
were obtained with the help of one-dimensional 
hybrid code CONEX [5, 6] in regime of cas- 
cade equations solution for interaction models 
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QGSJET 01 [7], SIBYLL 2.1 [8], NEXUS 3.97 [9] 
and QGSJET-II-03 [10]. In this paper all the re- 
sults are given for parameterizations of PCR nu- 
clei spectra from [1] with high helium flux. Nu- 
clei with A > 4 were treated in the framework of 
the superposition model, high accuracy of this ap- 
proach was confirmed by our calculations. More 
details of this check along with complete descrip- 
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Figure 2. Contribution of primary protons with 
energies E pr i m to the muon differential spectrum 
at sea level for ^ = 1.29 TeV. 



Figure 3. Inclusive spectra p + A — > 7r ± + X and 
p + A — > if ± + X (scaled down by 10) for incident 
proton with energy 10 TeV. 



tion of the calculation procedure will be given 
elsewhere. 

Comparison of the calculated muon fluxes with 
the experimental data, presented in Figure 1, re- 
veal familiar picture of high energy muon deficit. 
The reasons of its appearance were considered in 
our previous papers [11-13] and they still hold 
true regardless of the fact, that three more in- 
teraction models were included in our analysis. 
All interaction codes, except QGSJET 01, satis- 
factory describe data on muon flux only up to 
En ^100 GeV and then one by one fail to do it. 
Accounting that such muon energies correspond 
to primary energies above 1 TeV, studied with 
balloon(satellite)-borne emulsion chambers, one 
should simply relate muon deficit to underesti- 
mation of primary light nuclei fluxes, taking place 
in these experiments [11-13]. Unfortunately, dis- 
agreement between the models in the muon fluxes 
also appears at energies around 100 GeV, thus 
making impossible precise reconstruction of pri- 
mary nucleon spectrum for -E^rim > 1 TeV. In 
fact, in such conditions there are no reasons to 
rule out any of the models, except QGSJET 01, 
which, as it was said above, leads to remark- 
able disagreement with the experiment even in 
the range of reliable magnetic spectrometers data 



on PCR and muon spectra. 

To find why the models differ in the predicted 
muon fluxes let us consider quite characteristic 
energy of 1.29 TeV, where discrepancies between 
the models reach appreciable values and the data 
on muons from underground installations are yet 
quite reliable. Contributions of primary protons 
to the differential flux of muons of the given en- 
ergy, presented in Figure 2 show, that spread 
in muon fluxes between the interaction models 
is entirely due to uncertainties in the descrip- 
tion of 7T ± , _ftT ± -spectra in fragmentation region 
x = En^/Eprim > 0.1 (see Figure 3). Since in- 
clusive muon flux is sensitive nearly only to the 
characteristics of the very first primary particle 
interaction, hence, the harder these spectra are 
in the particular model, the larger muon inten- 
sity its use leads to. For the lower values of x, 
i.e. for i^prim > 10 TeV, all the models give prac- 
tically the same muon yields. As noted above, in 
view of uncertain situation with primary spectra 
for E pr i m > 1 TeV, one can not give preference 
to any of the models in comparison with the oth- 
ers. If to demand the minimal disagreement with 
the direct measurements data on PCR spectra, 
then obviously SIBYLL 2.1 satisfies this require- 
ment the best, or, in other words, one may say 
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Figure 4. Hadron spectra at the EAS-TOP depth 
t = 820 g/cm 2 . 



that it provides the most acceptable description 
of 7r ± , production spectra in p-air collisions 
in fragmentation region. 

3. Hadrons 

Comparison of our calculations with the most 
recent measurements of inclusive hadron flux, 
performed by EAS-TOP group [14], is presented 
in Figure 4. First, let us note the following 
facts. Below 100 GeV all the calculated spectra 
have breaks, because of non-perfect matching of 
low-energy model GHEISHA to the high-energy 
models. Shape of the measured hadron spec- 
tra also breaks at energies above 4 TeV and the 
data become less definite, thus in the forthcoming 
analysis we are going to use data only for ener- 
gies from 129 GcV to 4 TcV. For these energies 
QGSJET 01, QGSJET II and SIBYLL 2.1 quite 
reasonably reproduce the shape of the measured 
hadron spectrum, NEXUS 3.97 leads to spectrum 
with almost constant power index. One can see, 
that the most consistent description of the data 
for specified energies provide QGSJET 01 and 
SIBYLL 2.1. In contrast with the muons there 
are no energy range, where the models agree on 
the hadron fluxes and the reasons of this disagree- 
ment arc not as simply to point out as in the case 



Figure 5. Inclusive spectra tt^ + A — > n ± + X, 
p+A — > p+X (scaled down by 10), p+A — > n+X 
(scaled down by 50) for incident particles with 
energy 10 TeV . 



with muons. The most important characteristics 
in this analysis are total inelastic cross section, 
determining chances of primary particle to sur- 
vive, shapes of inclusive spectra p + A — > p + X, 
p + A -> n + X, 7T ± + A -> tt ± + X in the 
very forward region, responsible for substantial 
process of leading particles production (see Fig- 
ure 5 for the listed spectra) . Let us briefly outline 
the major conclusions, which one may come to in 
the given situation. NEXUS 3.97 gives the lowest 
fluxes as of hadrons in total, so of nucleons and 
mesons, and this happens in spite of the lowest 
inelastic cross-section values. Inclusive spectrum 
p + A — > p + X immediately helps to figure out, 
that incident protons in NEXUS 3.97 have com- 
parably low chances to save most of their energy 
in collision and this leads to such low nuclcon flux, 
the same may be said about meson flux and pro- 
duction of pions by pions. Similarly, from com- 
parison of the inclusive spectra, it can be easily 
understood, why QGSJET II gives the highest 
hadron flux. Note, that SIBYLL 2.1 concedes to 
QGSJET II in hadron intensity mostly because 
of less effective production of leading neutrons in 
p-air collisions and due to the larger total inter- 
action cross-section. 
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Thus, from analysis of the data on hadron flux 
it is difficult to imply any strict constraints on in- 
clusive spectra shapes, since mechanism of hadron 
spectrum formation is more sophisticated than 
that in the case of muon spectrum. SIBYLL 2.1 
and QGSJET 01 display quite a different be- 
haviours of the relevant inclusive spectra and to- 
tal interaction cross-sections, but both models al- 
most equally succeed in description of the EAS- 
TOP data. Alas, even this conclusion must be 
taken with care, since it is based on the single set 
of data and we have only indirect indications on 
the accuracy of this set, e.g. such as agreement 
of primary proton fluxes, obtained by E AS-TOP 
and KASCADE teams (the latter is derived from 
flux of unaccompanied hadrons [15]). 



muon fluxes for energies well below 100 GeV. 
Hence, one may conclude, that SIBYLL 2.1 gives 
the most acceptable overall description of the 
muon and hadron fluxes. Some underestima- 
tion of the muon flux for > 100 GeV almost 
for sure should be related to the underestima- 
tion of primary nucleon flux for E pT i m > 1 TeV. 
Overstatement of the hadron flux, which will be 
emphasized by this correction, may be compen- 
sated via slight reduction of diffractive events 
fraction. 

Authors greatly acknowledge CONEX team 
and personally Tanguy Pierog for their kind per- 
mission to use CONEX cross-section tables and 
for technical support. 



4. Self-consistency check and conclusions 

Self-consistency implies, that PCR mass com- 
position and spectra, once been retrieved from 
one kind of EAS data with particular interaction 
model, shall bring to satisfactory description of 
all other types of EAS data with the use of this 
very model. In our case this means, that one can 
reconstruct flux of primary protons, for example, 
from the data on hadron flux and then to apply 
it as input to get flux of muons. As we have seen, 
to match the data on muons with all the mod- 
els the primary nucleon flux for E pllm > 1 TeV 
must be increased in comparison with the direct 
measurements data. But, on the contrary, to 
describe EAS-TOP hadron flux with QGSJET II 
PCR flux must be significantly decreased, to less 
extent this applies also to SIBYLL 2.1. In the 
case of QGSJET II this change would result in 
very large disagreement with the data on muons. 
Since there are almost no need in correction 
of PCR spectra fits [1] for QGSJET 01 to agree 
with the data on hadron flux, this model also fails 
to satisfy self-consistency conditions. The only 
model, that simultaneously understates fluxes 
of hadrons and muons is NEXUS 3.97. But, 
even leaving aside misfit of the hadron spectrum 
shape, any changes in primary nucleon intensity 
for Sprim > 100 GeV/n, needed to minimize 
discrepancy with the EAS-TOP data, will im- 
mediately lead to disagreement with the data on 
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